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Oi 1 We discuss critical slowing-down of several gauge-fixing algorithms for the so-called 

fli ' A-gauges in the SU(2) case at zero temperature. For these gauges we also evaluate 

fT , ' the gluon propagator using different definitions of the lattice gluon field, corre- 

sponding to discretization errors of different orders. 

1 Gauge Fixing 

^ . EfRciency of gauge fixing algorithms is an important issue in the study of gauge- 

dependent quantities such as the gluon propagator on the lattice. We study 
^— V . critical slowing-down (CSD) by determining the dynamic critical exponent z 

^s) ' for various algorithms, using optimal gauge-fixing quantities. 

^^ . We consider the SU{2) case of the standard Wilson action and, in order to 

^^ ' impose the A-gauge conditionEI, we look for a local minimum of the functional 
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flj ' This corresponds to the condition 
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^ : £x[9] = -^Tr[Ag(x)t/4x)gt(x + erf) + ^ 5(x)[/^(x)5t(x -f e^)] . (1) 
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A5,4^)(x)+^a^4«)(x) =0, (2) 



H \ where A is the gluon field {U — t/^)/2i . The case A = 1 is the Landau gauge, 

whereas A = is the Coulomb gauge. At each step of a gauge-fixing algorithm, 
5(x) is updated to a new valufi.(7("'^"')(x) in order to minimize £[g\. We consider 
the following algorithms: crHu the local methods of Los Alamos (exact local 
minimization), Cornell (where the functional is minimized in the direction of 
the local downhill gradient of £[g\)^ overrelaxation (where local minimization is 
combined with an "energy-preserving" update by means of a tuning parameter 
w), and stochastic overrelaxation (where either type of update is randomly 
chosen at each step), as well as the global method of Fourier acceleration. This 
algorithm can be seen as a generalization of the Cornell method, whose update 
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is given by 

' 3(x) . (3) 



g("""')(x) ex i-a(V-A(9))(x) 
In fact, for the Fourier acceleration method we have 
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(new) / 



(x) (X h - «pL. A-1(V • A(«))(x)l g(x) . (4) 



The inverse of the A-Laplacian Aa = Xd^ + J2n=i ^t i^ usually obtained from 

l^^^p^^(k) F , where p\{}i) are the eingevalues of Aa, ancLF indicates the 
Fourier transform. In the Landau case we have introduced u a multigrid im- 
plementation of this method (where the inverse Laplacian is computed using 
a multigrid algorithm), showing improved behavior and applicability. 

In order to monitor the convergence of these algorithms, we evaluate at 
each step the quantity |V • A^^'p, as well as the quantity 

Q « EEE[Q^(^'^) -^^]'/[q:]' (5) 

U C Xi^ 

where Q^ ^ (1/7V.) E.„ Ql{^.) and Ql{x,) ^ E^^. E.,(^^^^)^(|) , which is 
seen to be particularly sensitive to the goodness of the gauge fixingS. For these 
two quantities we expect to observe an exponential decay with the number of 
iterations t. In the limit of large t we can introduce a relaxation time r such 
that the decay is given by exp( — t/r ) . As the lattice size N is increased, and 



at constant physics (see Section LI), r will grow as t ~ N^ , where z is the 
dynamic critical exponent of the algorithm. [Note that all the algorithms we 
consider, except for Los Alamos, need tuningJ 

In the Landau case we have simulated oa'u the above algorithms in the 2d 
case at finite values of the coupling /3, and in the 4d case at /3 = cxd. We have 
found that, as expected, the Los Alamos algorithm has z ^ 2, the improved lo- 
cal algorithms (i.e. Cornell, overrelaxation and stochastic overrelaxation) have 
z '^ 1, and tijat z ^ for Fourier acceleration. These results were also verified 
analyticallyQ at /? = oo. The multigrid implementation of Fourier acceleration 
is found to be considerably more stable and easier to tune than the standard 
method, and it can be used efficiently on parallel machines. Here we extend 
this study to the 4d case at finite (3 and to values of A different from 1 . 

1.1 Simulations and Results 

In order to analyze CSD for an algorithm, we have to evaluate r for different 
combinations of lattice size N and coupling /3, but at constant physics. This 
is done easily in d = 2 by keeping the product N^/a constant, since the string 



Table 1: Average number of iterations for the various gauge- fixing algorithms and for differ- 
ent A-gauges. We consider N = 8 at f3 = 2.6957. 
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FFTFA 


MGFA 


1 


631(42) 155(11) 137(12) 


159(8) 


95(11) 


94(10) 


0.5 


772(44) 176(21) 


143(9) 


178(10) 106(15) 


93(6) 



tension a is a known function of f3. No such function is available in d = 4. In 
this case we use the constant physics obtained H by keeping a lattice definition 
of the running coupling a^ constant. This prescription is well described by 
the fit /3 = 1.905 + 0.308 In A^ for /3 > 8. In particular, we consider A^ = 8 at 
/3 = 2.6957, N = 12 at P = 2.8485 and A^ = 16 at /3 = 2.9586. By fitting 
our data for the the relaxation time r to the function r = cN^ we obtain the 
dynamic critical exponents z fbijlhe various algorithms (neglecting finite-size 
effects). Our previous results BuQ are confirmed also in the four-dimensional 
case at finite values of /3: the Los Alamos method has z ~ 2, while the improved 
local algorithms have z ~ 1 and the Fourier acceleration method shows z ^ 0. 

All these results were obtained using the stopping criterion |V • Ap < 
10"^^ . We also check how the performance of the various algorithms is modi- 
fied if the condition Q < 10^^^ is used. For the case Ny^ 16 and /3 — 2.9586 we 
have found, in agreement with our previous findingsHiJ'cl, that for the stochas- 
tic overrelaxation and Fourier acceleration methods the number of sweeps is 
essentially independent of the quantity chosen for the stopping criterion, while 
the Los Alamos, Cornell and overrelaxation methods require about 20% more 
sweeps if the quantity Q is considered. 

Finally, in Table |l| we report the number of sweeps needed to satisfy the 
condition |V • A\'^ < 10~^^ , obtained for lattice size A^ = 8 and coupling 
P = 2.6957 for the various algorithms for two different values of A, namely 1 
and 0.5. From these data we see that the performance of the algorithms seems 
to depend weakly on the value of A. 



2 Gluon Propagator 

The study of the infrared behavior of the gluon propagator at zero temper- 
ature 19 provides a powerful tool for gaining insight into the physics of con- 
finement in non-Abelian gauge theories. In high-temperature QCD, the long- 
distance behavior of the gluon nrppagator is directly related to the electric 
and magnetic screening lengths Q'H. Here we study the gluon propagator at 
zero temperature, using different definitions of the lattice gluon field, and con- 
sidering different A-gauges. To this end, let us define the gluon propagator in 



momentum space as 
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D^ (fe) ^Y.[E cos(27rfci) Ql {t)]\ [^ sin(2^fct) Q^ (t)] ' (6) 



where Q^(t) = X^k « z ^/i(^' V^ -^' ^) ^'^'^ '^^ ^®* ^ ~ (*-*' '^j ^' ^)- -^^"^ ^'^^ lattice 
gluon field A^ we consider several possible definitions, leading to discretization 
errors of different orders. For example, we define: 

4..,,, , cyx)_^2w .„„ 4.,,., , i^«-'i' -^ Kc^^f ^ (7) 

If we set C/p(x) = exp[Jago <? • -^(x)], we obtain that both A); -^(x) and Ajj, '{x) 
are equal to ago <^ ■ -^(x) plus terms of order a'^ g^. We can also consider 

4=^)(x) EE 4/3AW(x) - l/3Af (x) (8) 

and 

AfCx) ^ 64/45 AW(x) - 20/454^)(x) + 1/45 I^^^iMl_Mf)^ . (9) 

It is easy to check that Al^ (x) = ago f? -^(x) plus terms of order a^ g^ , and that 
Af^ (x) = ago a ■ .A(x) plus terms of ordeiia^ g1. The definitions A^^^ and A'^'^^ 
were recently considered by Giusti et alEI. They find that the corresponding 
gluon propagators are equal modulo a constant factor. We perform a similar 
study here, evaluating D^^^\k) using the different definitions of the gluon field 
Aj^ given above, and we try to give an interpretation to this constant factor. 
Notice that Di' (k) + D2 {k) is related to the gluon propagator used by 
Karsch et alElB for the evaluation of the magnetic screening mass, and that this 
screening mass is invariant under rescaling of the propagators by a constant 
factor. 

2.1 Simulations and Results 

We have performed simulations at several values of j3, for lattice volumes V — 
8"*, 12"*, and for the gauge parameter A = 1 (Landau gauge) and A = 0.5. 
We obtain that, in all cases, the four propagators D'^*'(fc) are equal modulo 
a constant factor. Let us notice that this proportionality constant between 
different discretizations of the gluon propagator may be explained [at least for 
the simple cases D^^\k) and D(^)(fc)] as a tadpole renormalizationt2l. In fact. 



let us consider the tadpole-improved link C/^(x) = C/p(x)/uo, where uq is the 
mean link in Landau gauge. Then D^^' (k) gets multiplied by a factor Uq , and 
D'^^\k) by Uq^. At /? = 2.2 and A = 1 we have u^ = 0.6790(1). Thus, using 
tadpole-improved operators, the discrepancy £'(^^(fc)/Z?*^^)(fc) is reduced from 
1.862(9) to 1.264(6). (A similar analysis can be done at A = 0.5.) 

The gluon propagator is a gauge-dependent quantity. At finite tempera- 
ture, however, the pole masses obtained from the exponential decay of the 
gluon correlation functions at large spatial separations were proven to be 
gauge-independent, and therefore gluon propagators in different gauges should 
be related by a constant factor, as found in u. Here we consider only the 
zero-temperature case, and we expect to observe gauge dependence. We have 
checked that this is indeed the case, i.e. we observe gauge dependence for each 
given discretization D'^^\ Nevertheless, we see that, within a gauge, the be- 
havior of the gluon propagator is independent of the discretization modulo 
a renormalization factor. We plan to extend this study to the case of QCD 
at high temperature. A similar behavior in that case would mean that the 
screening masses are independent of the definition used for the lattice gluon 
field. 
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